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Abstract. The matrix equation XA + AX T = was recently introduced by De 
Teran and Dopico (Linear Algebra Appl. 434 (2011), 44-67) to study the dimen- 
sion of congruence orbits. They reduced the study of this equation to a number 
of special cases, several of which have not been explicitly solved. In this note we 
obtain an explicit, closed-form solution in the difficult Type 0-1 interaction case. 



1. Introduction 

The matrix equation 

XA + AX T = 0, (1) 
where A in M„(C) is fixed and X is unknown, was introduced in 2011 by F. De 
Teran and F.M. Dopico to study the dimensions of congruence orbits. In this 
setting, the codimension of the congruence orbit of A is given by the number of 
free parameters in the solution of (T). More recently the equation A} has also 
attracted the attention of a number of other authors @JJ6|. 

Now observe that for each fixed A, the solution space to (TJ is a Lie algebra, 
denoted g(A), equipped with Lie bracket [X, Y] — XY — YX. Building upon 
the initial work of De Teran and Dopico, a complete description of such twisted 
matrix Lie algebras is almost at hand. The computations in |3J provide explicit 
descriptions, up to similarity, of all possible twisted matrix Lie algebras except in 
a few highly problematic cases. In |6|, we provided an explicit description of the 
solution set in one of the most difficult cases. 

We are concerned here with the particularly troublesome setting where A is the 
direct sum of a so-called Type and Type I canonical matrix (see JZl Thm. 4.5.25], 
(3j Sect. 2]) These are matrices of the form /p(0) © T„, where / p (0) is the p x p 
Jordan matrix with eigenvalue and Tq is the q x q matrix 



r, ■, 



1 

-1 -1 

1 1 



(-1)1 







(-1)9+1- 

(-1)9 







2010 Mathematics Subject Classification. 15A24, 15A21. 

Key words and phrases. Congruence, congruence canonical form, Lie algebra, binomial coefficient, Pas- 
cal's Triangle, generating function. 

Partially supported by National Science Foundation Grant DMS-1001614. 

1 



2 



ALICE CHAN, LUIS ALBERTO GARCIA, STEPHAN RAMON GARCIA, AND AMY L. SHOEMAKER 



We aim here to provide an explicit, closed form solution to the matrix equation 

x(/p(o)er,) + (7 P (o)©r,)x T = o, (2) 

thereby obtaining a simple matricial description of the Lie algebra q(J p (0) © Tq). 
For even p, it is not difficult to show that 

8(7 P (o)er 9 )= fl (/p(o))e fl (r 9 ), 

where the Lie algebras fl(Jp(0)) and g(Tq) are elegantly described in |[3j p. 52, 
p. 54]. For odd p, however, the situation is much more complicated. Indeed, 
symbolic computation reveals that the matrices which arise are bewilderingly 
complex and typographically unwieldy. The consideration of this troublesome 
case takes up the bulk of this article. 

Besides providing a major step in the description of twisted matrix Lie alge- 
bras, our approach is notable for the use of a number of combinatorial identities. 
In fact, we conclude the paper by outlining an alternate description of the solu- 
tions to in which a variety of intriguing combinatorial quantities emerge. 



2. Solution 

We begin our consideration of the matrix equation ((2) by partitioning the 
unknown {p + cj) x (p + q) matrix X conformally with the decomposition of 
Jp(0) ®Tq. In other words, we write 

A B 
C D 



X 



(3) 



where A is p x p, B is p x q, C is q x p, and D is q x q. We should pause here to 
remark that the submatrix A just introduced is not the same as the fixed matrix A 
which appears in the general equation which plays no part in what follows. 

We next substitute the partitioned matrix into 10 to obtain the system 

A/p(0)+/ p (0)A T = 0, (4) 

DT q + TqD T = 0, (5) 

Br, + /p(0)C r = 0, (6) 

CJ p (0)+TqB T = 0, (7) 

of matrix equations which, taken together, are equivalent to 0. Fortunately, 
equations 10 and (0 are easily dispatched, since the explicit solutions in these 
cases have already been described by de Teran and Dopico 10 p. 52, p. 54]. 

It therefore remains to characterize the matrices B and C which satisfy © and 
10. This is by far the most difficult and involved portion of our work, although 
our task is made somewhat easier because B = C = when p is even |3, p. 60]. 

In the following, we assume that p is odd. Solving for B in 10 yields 

B = -} p (0) T C T T- T (8) 
and substituting this into © gives 

/p(0)C T = J p (Q) T C T Tq T Tq. (9) 
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Before proceeding any further, we require a few preliminary lemmas. The fol- 
lowing fact is well-known amongst those familiar with the congruence canonical 
form (see H p. 1016] or 1 9, p. 215]). 



Lemma 1. The cosquare Tq T Tq ofTq is given by 



where 





r, 


= ( 






"1 
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2 • • • 2 
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2 • • • 2 


A,:= 








1 ••• 2 










••• 1 



(10) 



(11) 



We also require the explicit reduction of Aq to its Jordan canonical form since 
the matrix which implements this similarity will be of great interest to us. Since 
we could not find the following result in the literature, we feel obliged to provide 
a detailed proof. 



Lemma 2. The matrix Aq satisfies Aq = Pq(lq + Jq(0))P« l , where Pq is the q x q 
upper triangular matrix whose entry is given by 



[P< 



2<t-i(-i)i-i 



Proof. Let Kq = Iq + Jq(0) and observe that i£j = l, then 

[A,P,] if i = S it 121- 1 = [P q Kq] iA , 

where 5u denotes the Kronecker .5-function, since the matrices involved are all 
upper-triangular. For j > 2, we make note of Identity 168 in [1], which states that 

whenever t > and n > 1. This gives us 

k=l 

k=i 

= [A 9 ]y[P 9 ]y+ E [A,]^^]^- 

k=i+l 

k=i+i \i ~ k ) 

= [P,]y+2^+ 1 E HO^f-"!) 

k=i+i V7 - V 
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[P,]y+2*-> +1 '£ (-1) 

/c=0 



k 



[P q } irj+2 i-i+\-iy 



(13) 
(14) 



= 

where the passage from l fT3l to l fl4|l follows by l(T2)> . Thus [A^P^; = [P^K^y for 
each pair (/,/), whence AgP^ = P^jR^, as claimed. □ 

We next require an explicit description of the inverse of Pq. For the sake of 
illustration, for £7 = 5 we obtain the corresponding matrices 



P 5 



~ 16 


-8 


4 


-2 


1 " 





8 


-8 


6 


-4 








4 


-6 


6 











2 


-4 














1 



1 

16 



11111 
2 4 6 8 
4 12 24 
8 32 
16 



Lemma 3. The inverse P q 1 of the matrix Pq is the q x q upper triangular matrix whose 
(i,j) entry is given by 



I- 1 



(15) 



Proof. Letting P' q denote the q x q upper triangular matrix whose entries are given 
by GD, we find that [P q P' q ]i,j = for ;' < i since the matrices involved are both 
upper- triangular. For j > i, we appeal to Identity 169 in (TJ, which we rewrite as 

" , / P 



k=0 



111 



"m,n 



for n > m, to conclude that 



k=l 



E 2 ^ fc (-i) 



k-i 



E(-ir 

it=i 



it-/ 



ifc-i 

k-i 



J 



7-1' 



— z/ \j — A: 



\Jfc— i 



EC- 1 ) 

fc=l 

£ ( _l)*-m 



k — i ) \j — k 



k=0 



k-i + lj \j -k 
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k=0 



k 

i - 1 



7-1 
k 



Thus P fl is invertible and P„ 1 = P'. 



Returning to © and applying Lemma [2] we obtain 

lp(0)C T P q 

In other words, the pxij matrix 
satisfies the equation 



J P (0)C T P q = (-i^jpiofcTpfy. 



Y = C 1 P a 



/ p (o)y = (-ir 1 / p (o) T y(z ( , + / ( ,(o)). 

Therefore the entries of Y satisfy 





^2,2 ■ 


■ yiA 


y3,i 


^3,2 ■ 


■ VM 




Vp,2 ■ 


■ VVA 












-1) 



q-1 





yu 

J/2,1 



yi,i + yi,2 
y2,i + y2,2 



.yp-i,i y?-i,i + yp-1,2 



o 

yi,ij-i + i/i,<? 

V2,q-\ +yi,q 
Vp-lj-l + yp-l,<;. 



from which we see that 



yy = HO' 1 {yi-i l j-i+yi-2,j) 



holds for 3 < i < p and 2 < j < q. This prompts the following lemma. 
Lemma 4. Ifi is even, then = Ofor 1 < j < q. If i is odd, then 



-(<?-!) 



□ 



(16) 
(17) 

(18) 

(19) 



(20) 



Proof. Comparing entries in il8l , we find that 1/2,; = for 1 < j < q. In light of 
(Hill, it follows that 

y 4j - = (-l)^ 1 (y 2 , ; -i+y 2/ ;) =0, 
for 1 < j < q. Proceeding inductively, we see that y, y = whenever i is even. 

Now suppose that i is odd. For the basis of our induction, observe that d20l 
holds trivially when i + j = 2 (i.e., when i — j — 1). Now suppose that n > 3 and 
that | |2"0l holds if i + j = n — 1. Under this hypothesis, we wish to show that d20l 
also holds if z + / = n. In light of | fl9t , we have 

yy = (-i)' ?_1 (y ! --2,;-i+yi-2,;') 



-1) 



'/-l / i-3 

50-T-i 



; / t-3 



yu 
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-!)¥(<?-!) 



E 

k=l 

i 



i-3 
2 

-k-1 



i-3 
2 



yi,k 



= (-i)¥(«?-D £ 

fc=l v 

which completes the proof. 



i-1 

2 



2/l,fc> 



□ 



Armed with an explicit, entry-by-entry description of the matrix Y = C T Pq, 
we can construct a convenient basis for the solution space of fiTfr . Noting that 
Y depends upon the ^ free parameters y\ r \, J/i,2/- ■ ■ we introduce a family 
Wj, W2, . . . , Wn of matrices by letting W<j denote the matrix Y which corresponds 
to the initial values yij = Sj^ along the first row and whose remaining entries are 
determined by EB . To be more specific, it follows from < f20t that 



i-l 



if i is even, 
if i is odd. 



(21) 



In general, for i odd we have 

mh = 



l)¥(?-i) 



i-l 
2 



[W- 



iJy-(^-i) 



Wi(Mo))'" 1 



i-l 
2 



so that 

For example, if p 

Y = 
and 



*V< = Wi(/?(0)) 
5 we obtain 



£-1 



(22) 
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2yi,i + yi,2 


3/1,1 + 2yi,2 + yi,3 


yi,2 + 21/1,3 + 3/1,4 


yi,3 + 2yi, 4 + yi,5. 
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1. 



By ((21) we can see that [W,]^ = <5/ y for 1 < i,j < q whence the matrices 
Wi/ W2, . . . , are linearly independent. Considering j22)l . we now set W = Wi 
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so that for any solution Y of fl7t we have 



for some constants ax, «2/ . . . , a,. 

Returning to fhSl l and solving for C yields 



£-1 



(23) 




(24) 



where 1 is given explicitly by flHI l and W = Wi is given by J2T) . Turning our 
attention now toward ((8), we find that 



B = -/ p (o) T c T r- 



(25) 



so that B is completely determined by C. In particular, the dimension of the 
solution space of 10 and 10 is q, which is consistent with |3, Lem. 10]. Further- 
more, |3, Lem. 4] and |3j Lem. 5] tell us that the solution sets of ((4j) and 10 have 
dimensions [£] and respectively. Putting everything together, we conclude 
that the dimension of the solution space of our original equation $2$ for odd p is 



q + 



J- 



Appendix A. Combinatorial interpretation 



Although we are now in possession of an explicit description of the solution 
to l|2), it is worth examining the submatrices B and C a little bit closer since the 
entries of these matrices display some remarkable combinatorial properties. 

Using our closed form solution |24}, let us continue with our previous series 
of examples corresponding to the parameters p = q = 5. In this setting, we have 



C 



H 
16 

K(«l+4(«2 + «3)) 
^ («i +6a 2 +12«3 + 8*4) 



3- 

u 16 

1 L(3« 1 +2a 2 ) 

^(5*! +8*2+4*3) 

^ (7*i + 18x 2 + 20n 3 + 8a 4 ) 

^. + 2a 2 + ^a+3a 4 +a 5 



^(5*i+2«2) 
^ (13«, +4(3« 2 + « 3 )) 

^ (25a! + 38a 2 + 28ff 3 + 8a 4 ) 
41a n , ll!t 7 , 13flo , . 



In particular, the simple parametrization at the heart of l[24} is obfuscated by pre- 
and post- multiplication with the somewhat complicated matrices P^ T and W T , 
respectively. This results in the rather overwhelming complexity apparent in the 
preceding example. However, upon performing an appropriate reparametriza- 
tion, we can obtain the significantly simpler matrix 



C = 



Pi 
ft 
ft 

ft 
ft 



o ft o 

ft + 2ft 

ft + 2ft + 2ft 

ft + 2ft+2ft+2ft 

ft +2ft + 2ft+2ft +2ft 



ft 

ft + 4ft 
ft + 8ft + 4ft 
ft + 12ft + 8ft + 4ft 
3 5 + 16ft + 12ft + 8ft + 4ft. 



(26) 



s 
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We sketch here an independent combinatorial argument for determining the exact 
coefficients which arise in this manner. Since our motivation for doing so is purely 
aesthetic, we leave many of the tedious details to the reader. 
We begin anew, expanding |(9]l and using Lemma Q] to obtain 



Cl,2 
c l,3 



Cl,p 





C2,2 
c 23 



C2,p 





C?,2 
S3 

c q,p 




(_l)f-l 





ci,i c 2 ,i + 2c w 

^1,2 C 2/2 + 2c 1/2 

_ c l,p-l c 2,p-l + 2ci, p _i 







At this point, a short inductive argument tells us that c, y = whenever y is even. 
On the other hand, the odd columns of C are governed by a much more intricate 
recurrence. We begin by noting that there are no restrictions placed upon the 
entries in the first column of C. We therefore set 

c\,\ = j6i, c 2 ,i = j6 2 , . . . , c ?; i = j6 ? , 

where f>\, [$2, ■ ■ ■ , fiq are free parameters. The remaining odd columns can be 
computed using | |27|| , which provides us with the recursion 

aa+i = + E 2c ^0- (28) 

From this one can deduce that the entries of C satisfy 




where Sjt(n) denotes the coordination sequence for a fc-dimensional cubic lattice 
(see Table [TJ. To be more specific, Sjt(w) represents the number of lattice points in 
at distance n from the origin. For instance, the first row of Table [1] is 1,0, 0, . . . 
since the trivial lattice contains only the origin. The second row of Table [T] is 
1,2, 2, . . . since the lattice Z 1 = Z contains exactly one point at distance one from 
(namely itself), and exactly two points at distance n from (namely in) for 
n > 1. See [2 J for more information about the numbers S; c (n). 

Although we do not indulge in the somewhat lengthy derivation of d29t , it is 
possible to highlight the main points. First observe the Toeplitz matrix A t; , given 
by fTT) . which appears in |9]l is a truncation of the infinite Toeplitz matrix having 
symbol 

1 + X <•) n 

= 1 + 2x + 2x + 2x H . 

1 — x 

This gives rise to the coefficients evident in the recursion | |28l . More generally, we 
have the generating function identity 

(\^) k =t S k(n)x\ (30) 

which can be used to explain the appearance of the S^(n) which arise through 
repeated applications of the recursion (|28l . 
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n 
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n 
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n 
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n 
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Si(n) 


1 


2 


2 


2 


2 


2 


2 


2 


2 


2 


S 2 (n) 


1 


4 


8 


12 


16 


20 


24 


28 


32 


36 


S?,(n) 


1 


6 


18 


38 


66 


102 


146 


198 


258 


326 


Si{n) 


1 


8 


32 


88 


192 


360 


608 


952 


1408 


1992 


S 5 (n) 


1 


10 


50 


170 


450 


1002 


1970 


3530 


5890 


9290 


S 6 (n) 


1 


12 


72 


292 


912 


2364 


5336 


10836 


20256 


35436 


S 7 (n) 


1 


14 


98 


462 


1666 


4942 


12642 


28814 


59906 


115598 


Ss{n) 


1 


16 


128 


688 


2816 


9424 


27008 


68464 


157184 


332688 


Sg(n) 


1 


18 


162 


978 


4482 


16722 


53154 


148626 


374274 


864146 


S w (n) 


1 


20 


200 


1340 


6800 


28004 


97880 


299660 


822560 


2060980 



Table 1. The initial terms of the sequences So(w), Si (n), . . . , Sio(w). 
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G (n) 
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1 


1 


1 


1 


1 


d(n) 


1 


3 


5 


7 


9 


11 


13 


15 


17 


19 


G 2 (n) 


1 


5 


13 


25 


41 


61 


85 


113 


145 


181 


G 3 0) 


1 


7 


25 


63 


129 


231 


377 


575 


833 


1159 


G 4 (n) 


1 


9 


41 


129 


321 


681 


1289 


2241 


3649 


5641 


G 5 (n) 


1 


11 


61 


231 


681 


1683 


3653 


7183 


13073 


22363 




1 


13 


85 


377 


1289 


3653 


8989 


19825 


40081 


75517 


G 7 (n) 


1 


15 


113 


575 


2241 


7183 


19825 


48639 


108545 


224143 


G s (n) 


1 


17 


145 


833 


3649 


13073 


40081 


108545 


265729 


598417 


Gg(n) 


1 


19 


181 


1159 


5641 


22363 


75517 


224143 


598417 


1462563 


Gio(rc) 


1 


21 


221 


1561 


8361 


36365 


134245 


433905 


1256465 


3317445 



Table 2. The initial terms of the sequences Gq(«), G\(n), . . . , G\o(n). 



Recalling that B = -} p (0) T C T r- T one can rapidly deduce that the odd in- 
dexed rows of B are identically zero. For the even indexed rows of B, we may 
appeal to the generating function identity 

■ T^- x = ( t SkW A (LA = t Gk(n)*», 

1 */ 1 x \n=0 / \n=0 / «=0 

where Gfc(n) = E" = o %(^) denotes the number of lattice points in Z ,c at distance 
< j? from the origin (see Table O, to eventually obtain 

{0 if i is odd, 

(_l)/+(f+D(^) £ G.^W^, i£i = 2ik. 
i=0 
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For example when p = 7 and q = 5, we have 



B 



o 

-(05 + 04 + 03 + 02 + 0l) 

o 

- (05 +304+ 5ft;+7fe+9fr) 


(05+504 + 1303+2502+410!) 




04 + 03 + 02 +01 


04 +303 +502 +70! 


04 + 503 +1302 +2501 








(03 + 502+130i) 




(03 +302 +50i) 





(03 + 02 + 01) 













-i 



. 
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